Non-Fraunhofer Interference Pattern in Inhomogeneous Ferromagnetic Josephson Junctions 
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Generic conditions are established for producing a non-Fraunhofer response of the critical supercurrent sub- 
ject to an external magnetic field in ferromagnetic Josephson junctions. Employing the quasiclassical Keldysh- 
Usadel method, we demonstrate theoretically that an inhomogeneity in the magnitude of the energy scales in 
the system, including Thouless energy, exchange field and temperature gradient normal to the transport direc- 
tion, influences drastically the standard Fraunhofer pattern. The exotic non-Fraunhofer response, similar to that 
observed in recent experiments, is described in terms of an intricate interplay between multiple '0-7r'-states and 
is related to the appearance of proximity vortices. 



PACS numbers: 74.50.+r, 74.45,+c, 74.78.FK 

The well-known Fraunhofer diffraction pattern of the 
critical Josephson current has been extensively studied in 
superconductor/normal-metal/superconductor (S/N/S) junc- 
tions [1, 2]. The interest in how a supercurrent responds to 
an applied magnetic flux derives from the fact that this prop- 
erty is the key element in ultra-sensitive devices such as su- 
perconducting quantum interference devices (SQUID) [3-5]. 
Whereas S/N/S junctions are known to display Fraunhofer 
diffraction in the wide junction limit, the critical current de- 
cays monotonically as a function of the applied flux in the 
narrow junction limit. The crossover between these two dis- 
tinct types of behavior was theoretically described in terms of 
proximity vortices in the normal wire [6]. 

More recently, the orbital response of the supercurrent 
in magnetic Josephson junctions has attracted much interest 
[3, 4]. When the normal-metal interlayer is exchanged with 
a ferromagnet, thus forming an S/F/S junction, a new mech- 
anism comes into play compared to the S/N/S case. The 
ground-state phase difference between the superconducting 
reservoirs may then take the value of or it, depending on pa- 
rameters such as temperature and ferromagnetic barrier thick- 
ness [7]. Not only does this cause the supercurrent in magnetic 
Josephson junctions to decay in a non-monotonic fashion, but 
it was recently reported that non-Fraunhofer interference pat- 
terns appear in S/F/S junctions composed of alternating 0- and 
7r-states [3, 4]. Whereas the supercurrent is maximal at zero 
flux in the non-magnetic case, the supercurrent instead dis- 
played a minimum at zero flux in the S/F/S case [4]. These ex- 
perimental findings have motivated theoretical investigations 
[5, 8]. Non-Fraunhofer interference patterns have also been 
studied in S/I/S junctions with arrays of resistors [9]. 

Motivated by this, the following question is answered in 
this Letter: under which general conditions may the criti- 
cal supercurrent respond to an external magnetic field in an 
anomalous fashion, producing a non-Fraunhofer interference 
pattern? We establish these conditions and moreover explain 
the origin of this exotic phenomenon. To do so, we solve the 
quasiclassical Keldysh-Usadel equations. In the majority of 
past works, the investigation of the non-Fraunhofer patterns 
were restricted to incorporating a linear ansatz for the be- 




FIG. 1: Experimental setup of the inhomogeneous ferromagnetic 
Josephson junction. The green arrow shows the direction of the 
inhomogeneity in magnitude of either the Thouless energy, mag- 
netic exchange field or temperature normal to transport direction (x- 
direction). To model an inhomogeneity in the Thouless energy, the 
setup A) is considered in this letter. Although making the set up B) 
might be easier to achieve technically, the two setups generate the 
same results in the diffusive limit. An inhomogeneity in the mag- 
netic exchange field is modeled by h = h(0, 0, y/W). The external 
field is applied to the system in the z-direction (not shown). 



havior of the superconducting U(l) phase [3, 4]. In contrast, 
we have employed in this Letter a model of a ferromagnetic 
Josephson junction which takes into account an external mag- 
netic field with an arbitrary dependence on the coordinates and 
direction of the field. This model allows us to study the critical 
supercurrent through an inhomogeneous junction without re- 
course to any ansatz. The possibility of having an arbitrary in- 
homogeneous magnetization texture in the F region makes the 
model highly general. The results of the developed theory are 
qualitatively in agreement with the recent experimentally ob- 
served non-Fraunhofer patterns. Remarkably, we find that the 
critical Josephson current through the F region is suppressed 
at zero-external magnetic field within the wide junction limit 
when the magnitude of any of the energy scales of the system, 
i.e. Thouless energy, exchange field and temperature are in- 
homogeneous normal to the transport direction. Crucially, to 
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achieve a non-Fraunhofer response, the inhomogeneity must 
include at least one 0-7r-state. In this case, the second peak of 
magnetic interference pattern of critical supercurrent becomes 
larger than the first, in contrast to the Fraunhofer pattern. We 
explain this behavior in terms of 0-tt crossover states and also 
relate our results to the appearance of proximity vortices in- 
side the F region. 

Consider the schematic of the proposed experimental setup 
in Fig. 1. The inherent Josephson penetration depth Aj is as- 
sumed to be larger than the width of junction, such that one 
may avoid screening effects imposed by the Josephson cur- 
rent on the external magnetic field. This field is assumed to 
be directed along the z-direction. We work with a vector po- 
tential satisfying the Lorentz gauge i.e. \^-A=0 and choose 
specifically A=—Hyx in which H represents the strength of 
external magnetic field. The magnetic flux due to the intrin- 
sic magnetization of the ferromagnetic region is ignored, as 
is known to be a good approximation in most cases [4]. To 
investigate the transport properties of this system, the quasi- 
classical theory of superconductivity in the diffusive regime 
is employed, so that the Gor'kov equations are reduced to the 
Usadel equations [10]. The Usadel equation inside the F re- 
gion together with appropriate boundary conditions is used for 
obtaining observable quantities of the system. In the presence 
of a static external magnetic field, the Usadel equation is suc- 
cinctly given by; 

D[d,g(x,y,z)[d,g(x,y,z)]} + i[ep 3 + 

diag[h(x, y,z) -a, (h(x,y,z) ■ a) T },g(x, y, z)\ = 0,(1) 

where h(x, y, z) stands for exchange energy, g is the full 8x8 
Green's function matrix, while p 3 and a are 4 x 4 and 2x2 
Pauli matrixes, respectively [18]. Here D is the diffusion con- 
stant and d = VI — ieA(x, y 7 z)p 3 . Within the weak prox- 
imity regime, one may can expand the Green function around 
the bulk solution g i.e. g(x, y, z) ~ g 4- f{x, y, z), where 
g = diag(4,-l_) [11]. Therefore, the retarded component of 
Green's function reads: 
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The advanced and Keldysh blocks are also 
given via g A (x, y, z)=-(p 3 g R (x, y, z)p 3 )* and 
g K (x, y, z)=(g R (x, y, z)-g A {x, y, z))tanh(e/2fc B T) un- 
der equilibrium conditions. At the two N/S interfaces the 
Kupriyanov-Lukichev boundary conditions [12] are com- 
pactly written by: 2Cg[(V* - ieAp 3 ) ■ n,g] = [g BCS (0) , g], 
the ratio between the resistance of the barrier region and the 
resistance in the F film is defined as (=Rb/Rf, gBCs(^) * s 
the Green's function in the two superconductor reservoirs 
[13] and n is a unit vector normal to the interface [11]. At the 
vacuum borders, the Green's function satisfies d y g = 0. 

The Usadel equations in their present form constitue a set 
of complicated coupled differential equations which we have 
solved numerically by using a collocation method. Thus, 



the approximate solution components are assumed to be lin- 
ear combinations of bicubic (tricubic, for three-dimensional 
problems) Hermite basis functions, and required to satisfy the 
Usadel equation exactly at 4 (8, for three-dimensional prob- 
lems) collocation points in each subrectangle of a grid, and 
to satisfy the boundary conditions exactly at certain bound- 
ary collocation points [14]. Finally, Newton's method is 
used to solve the (nonlinear, generally) algebraic equations 
resulting from the collocation method formulation [17]. In 
order to study transport properties of the inhomogeneous 
junction, the current density through the junction is con- 
sidered: J(i?,0) = J JdeTi{p 3 (g{d,g]) K }, here J = 
N^eD/A and jV is the number of states in the Fermi sur- 
face. Performing an integration over the y-coordinate pro- 
vides the total supercurrent flowing through the junction, 
I(4>) = Iq J J dydeTr{p 3 (g[d,g]) K }. To understand the 
magnetic interference patterns of such junctions, we also in- 
vestigate the spatial variation of pair potential inside the F re- 
gion calculated via: U — [7oTr{(/5i — ip^) J der 3 g K }, where 
Uq = —N \/16 [18]. The temperature, width and lower base 
of the wedged junction are fixed at T/T c = 0.01, W/£s = 10 
and d/(_s = % respectively (the so-called wide junction limit). 
The proximity controlling parameter ( is also fixed at 5, en- 
suring that we operate in the weak proximity regime. Energy 
units are used so that h—kB=\. 

The results for the critical Josephson current through the in- 
homogeneous S/F/S junction as a function of normalized ex- 
ternal magnetic flux i.e. are presented in Fig. 2. In 
frame A), the magnitude of the Thouless energy is inhomo- 
geneous in the y-direction: the F region has a wedged shape 
[see A) in Fig. 1]. In frame B), the magnitude of magnetic ex- 
change interaction is inhomogeneous in the y-directions and 
follows a h = h(0, 0, y/W) pattern. The normalized critical 
current IJIq exhibits a suppression at zero external flux for 
some values of the wedge angle a. In the case of trapezoidal 
junction, the second peak in the interference pattern takes a 
larger value than the first for an interval of a-values. Typ- 
ically, this behavior is enhanced at a=7r/43 and then disap- 
pears for larger values a > ir/20 (see inset panel of panel A) 
of Fig. 2). A similar magnetic interference pattern is gener- 
ated when the magnitude of exchange field is inhomogeneous, 
as shown in frame B). In this case, the non-Fraunhofer pat- 
tern phenomenon is pronounced for e.g. h = 4. 27 A. The 
results show qualitatively good consistency with recently re- 
ported non-Fraunhofer patterns for "0-7r"-stacks in Ref. 4. For 
reasons to be described below, we expect that the same non- 
Fraunhofer magnetic pattern would arise when the tempera- 
ture of the system along the y-direction is variable and has an 
inhomogeneous form. We have also investigated (not shown) 
other magnetization textures such as domain-wall, skyrmion 
and spiral (with helical axis normal to transport direction), 
and found that they generate the standard Fraunhofer patterns 
because of the constant magnitude of the magnetic exchange 
field, |h|=h. 

In comparison, the behavior of critical Josephson current 
through a S/N/S wedged junction is investigated in Fig. 2 C). 
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FIG. 2: Normalized supercurrent through an inhomogeneous Josephson junction vs normalized external magnetic flux $/$o perpendicular 
to the junction. A): An S/F/S junction where the F region has a wedged shape with inclination angle a and exchange field h = 10 Ao. The 
inset panel zooms in on the interference pattern at small and large values of $ and a, respectively. B): The F region is now geometrically 
rectangular, but the magnitude of exchange field is now inhomogeneous according to the texture h = h(0, 0, y/W). C): An S/N/S junction 
where the N region has a wedged shape with inclination angle a. The legends are the same as in A). In all cases, the height of the trapezoidal 
region is fixed at W = 10£s while the upper base d is equal to 2£s (see Fig. 1). Therefore, a=0 makes a rectangular junction (L — d — 2£s) 
and tt/130 (L = 2.483£ s ), tt/67 (L = 2.938£ s ), tt/43 (L = 3.463£s) and tt/20 (L = 5.168£s) make wedged junctions. 



For a = 0, we recover the results of Ref. 6. With increasing 
a, the normalized supercurrent is subject to an overall reduc- 
tion, because the effective junction length increases due to the 
inhomogeneity in the magnitude of Thouless energy. Unlike 
the inhomogeneous S/F/S case above, however, the first peak 
in the diffraction pattern is larger than others for all values of 
a, exhibiting the standard Fraunhofer pattern. Therefore, the 
exotic non-Fraunhofer pattern only can be observed in ferro- 
magnetic junctions under the conditions discussed above. 

To further understand the outstanding difference between 
interference patterns of the homogeneous and inhomogeneous 
ferromagnetic junctions, we consider how the presence of an 
external magnetic field influences both the 0-ir transition pro- 
file of the S/F/S junction and the proximity vortices pattern 
in the F region. The origin of the suppressed central peak in 
the S/F/S wedged Josephson junction (a ^ 0) is mainly stud- 
ied in this Letter, and we argue why this mechanism accounts 
for the non-Fraunhofer pattern in other cases where there is a 
magnitude gradient of the exchange field and/or temperature 
along the direction normal to the transport direction so that it 
includes at least one 0-7r-state. To this end, we will later in- 
vestigate the current density spatial map of magnetic junctions 
and compare a=0 with a / 0. Fig. 3 reveals an illustrative 
profile of the normal and ferromagnetic Josephson junctions. 
Part A) illustrates a spatial map of pair potential in the normal 
region of the S/N/S junction where a=0 and tt/43 for the left 
and right panels, respectively. The upper base of the trape- 
zoidal region is fixed at d — 2£s while for a=7r/43 the lower 
base takes the value i=3. 463£s. The increment of a deforms 
the proximity vortex pattern compared to the pattern of the 
rectangular junction. The distance between two neighboring 
vortices is no longer equal to Q$H/ d in contrast to that of rect- 
angular junction within the wide junction regime (left panel). 
The spatial maps of the pair potential are given for $=4$o 



and zero superconducting phase difference i.e. (f>=0. Now, the 
increment of a removes the proximity vortices inside the nor- 
mal segment of the junction gradually. The variation of super- 
conducting phase difference tj> however, moves the vortices 
along the y-direction in both the normal and ferromagnetic 
junctions [6]. Part B) exhibits the equivalent investigation for 
a ferromagnetic Josephson junction. In the rectangular case, 
the pair potential shows the same behavior as the vortex pat- 
tern as S/N/S case. In contrast, the vortex pattern is highly 
deformed in the wedged ferromagnetic junction. A zoom-in 
is shown for the middle of F wire with a=7r/43 using differ- 
ent color map. The strong deformation may be understood by 
noting that the increment of a effectively synthesizes multi- 
ple "0-7r" states in the same junction and also from the pair- 
breaking of the exchange field. When a becomes non-zero, 
the junction may be thought of as a superposition of multiple 
and 7r junctions. 

To understand this quantitatively, parts C) and D) should 
be considered together. Part C) illustrates the Q-tt crossover 
profile where a = for two different values of external flux 
$=0, 3$o/2 as a function of F-layer length d/£s- The first 
and second transitions occur at <i=2.38£,s and 3.46£s, respec- 
tively. It is important to note that the latter length is identical 
to the lower base of the trapezoidal junction when a=7r/43. 
The plot also shows that applying an external magnetic field 
reduces the magnitude of the current nonlinearly, although 
the locations of 0-tt points are left unchanged. The anoma- 
lous Fraunhofer diffraction pattern of the critical supercurrent 
can now be well understood by noting part D). Increasing the 
junction angle a renders more parts of the junction to have op- 
posite supercurrent flow direction which then partially cancel 
each other. One should note that in the trapezoidal region, the 
amplitude of critical current is non-uniform. More ir -parts, 
therefore, are needed to cancel the 0-parts of the junction that 
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FIG. 3: A) Normalized spatial map of the pair potential for a S/N/S junction where a=0, 7r/43 (L=3.463£s) in the left and right frames, 
respectively. B) The normalized pair potential in a S/F/S junction. Arrows indicate the location of proximity vortices. The zoom-in frame of a 
S/F7S trapezoidal junction with q=7t/43 is shown using a different color map. The magnetic flux through the N or F region is assumed to be 
<I>=4$ and no superconducting phase difference is applied (4>=Q). C) The 0-7T crossover profile of a rectangular S/F/S junction where $=0, 
3<I?o/2 and W = 10£s vs. normalized junction length d/£s. D) Current density spatial map of the magnetic Josephson junction in the absence 
of external magnetic field, $=0. Top and bottom frames exhibit current density flowing through the junction where a=0, 7r/43, respectively 
(a=7r/43 constitutes a '0-7r'-junction). Blue arrows indicate current directions. 



occur for the top region with smallest effective length L. A 
key observation is that the above results suggest a venue for 
producing an anomalous non-Fraunhofer interference pattern 
resulting without necessarily distorting the geometry of the 
system, in effect allowing for anomalous interference even in 
standard rectangular junctions. 

In conclusion, we have proposed experimentally accessible 
generic conditions for achieving non-Fraunhofer interference 
patterns of the critical supercurrent as a function of external 
magnetic flux. The key property is the controllable numbers 
of gradual 0-ir states in the same junction by incorporating 
an inhomogeneity in the magnitude of energy scales of sys- 
tem i.e. Thouless energy, exchange field and/or temperature 
normal to the transport direction. We examine the proposed 
generic conditions for some limiting cases and find good qual- 
itative consistency with the recently observed non-Fraunhofer 
magnetic interference patterns in 0-ir stacks. 
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